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ABSTRACT 

Binaries are very common in galaxies, and more than half of Galactic hot sub- 
dwarf stars, which are thought as a possible origin of UV-upturn of old stellar 
populations, are found in binaries. Previous works showed that binary evolu- 
tion can make the spectra of binary star populations significantly different from 
those of single star populations. However, the effect of binary evolution has not 
been taken into account in most works of spectral fitting of galaxies. This paper 
studies the role of binary evolution in spectral fitting of early-type galaxies, via 
a stellar population synthesis model including both single and binary star pop- 
ulations. Spectra from ultraviolet to optical band are fitted to determine a few 
parameters of galaxies. The results show that the inclusion of binaries in stel- 
lar population models may lead to obvious change in the determination of some 
galaxy parameters and therefore it is potentially important for spectral studies. 
In particular, the ages of young components of composite stellar populations be- 
come much older when using binary star population models instead of single star 
population models. This implies that binary star population models will measure 
significantly different star formation histories (SFHs) for galaxies compared to 
single star population models. In addition, stellar population models with binary 
interactions measure larger dust extinctions than single star population models 
on average. It suggests that when using binary star population models instead 
of single star population models, negative extinctions are possibly unnecessary 
in spectral fitting of early-type galaxies. 



Subject headings: Galaxies: elliptical and lenticular, cD — Galaxies: stellar 
content — Galaxies: star formation. 
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Introduction 



Binaries are common in galaxies , and their evol ution is often different from single stars 
if binary components are not too far (jEggletorul2006[ ). In our Galaxy, the fraction of binaries 
is shown as large as about 50%. Binaries are found in different parts of the Galaxy, e.g.. 
Gala c tic field, halo, and thick and thin discs (Ipuauennov fc Mavorl llQQlt iFischer &: Marcy 
1992 : Mayor et al. 1992 : Raghavan et al. 201ol: Carney et al. 20051 ). In addition, most stars 
massive than the sun are born in binaries (IBrinchmannl l2010l ). and the massive the stars, 
the larger the fraction of binaries ( iLada fc Ladall2003l ). Therefore, binaries are absolutely 
important for understanding the stellar populations of galaxies, and it is necessary to study 
the role of binary evolution in stellar population studies. The importance of binaries has 
been recently lis ted as one of six irn portant challenges in stellar population studies in the 
next decades by IBrinchmannl ( 120101 ). However, most previous works used single stars for 
modeling stellar populations because of the complexity of binary evolution and more input 
parameters in binary star stellar population synthesis (bsSPS). Besides, astronomers hope 
that binary evolution does not change the results of single star stellar population (ssSP) 
studies too much. If this is true, to take single star stellar population synthesis (ssSPS) will 
be a good choice for stellar population studies. 

i et all bood : iHan et all boO?! : \u fc Hani boOSal M 



Meanwhile , a few works (e.g.. 



Li et al.l l2012at [Hernandez fc Bruzuall 1201 ll ) have tried to investigate the effects of binary 



evolution on stellar populati on studies and have used bi n ary star stellar populat i on (b sSP) 
model in a few works (e.g.. 



Li et al. 


2006. 


2010. 


2012b: 



shows that binary evolution is able to reprod uce many interesting results such as blue strag- 



glers, red stragglers, extended main sequence Jhi fc HanlboOSal iLi et al 



excess (hereaft er UV-upturn) in elliptical galaxies (e.g., iHan et al.l bOOTt iLi fc Hanll2008a 



Li et al. 



1994 



2012b), and UV flux 



2012a[). without any special assumptions like very old age or high mass-loss (see 



Lee 



Park &: Leelll997l : iBressan et al.lll994l : lYi et al.lll997l for comparison) . Binary evolution 



of Galactic 


lot sub dwarf stars are found in binaries ( 


Ferguson et a 


— 
. 1984 


: AUard et al. 


1994: 


Theill et al. 


1995: U] 


la & Theill 1998: 


Aznar Cuadrado &; Jefferv '. 


2001 


: ] 


Maxted et al. 


2001: 


Williams et al. 


2001; 


Reed & Stiening 


2004). and it is well established that the vast majority 



(and quite possibly all) of hot subdwarfs are the results of binary interactions (IHan et al. 



20071 ). On the other side, the widely used ssSP models fail to reproduce UV-upturn for old 
stellar populations under normal assumptions. Thus it suggests that binaries are important 
and possibly the main contributors to UV-upturn spectra of old populations. In addition, it 
shows that binaries can change the results of simple stellar population (SSP) studies obvi- 
ously (ILi fc Hanll2008bl ). when spectral line indices or colours are used to determine stellar 
population parameters (age and metallicity) of galaxies. Although the importance of bina- 
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ries has been shown, it is far from well understanding the role of binaries. Deeper and more 
detailed studies are therefore needed. Because it is shown tha t binary evolut ion leads to 
significant change in the UV and optical spectra of populations (ILi et al.ll2012af). and fitting 



such spectra are the bas is of most stellar populatio n stud i es (see, the rev i ew of 



2011 



or papers such aslCid Fernandes et al.l 12004 l2005l: lLu et al.ll2006t iTojeiro et al. 



Koleva et al.ll2008l : iTojeiro et al.ll2009l : iKoleva et al.ll2009l ). it is important to study the effect 



Walcher et al. 



2007 



of binaries on full spectrum fitting. 

This work aims to investigate the potential importance of binary evolution in spectral 
studies of early-type galaxies. The reason for studying early-type galaxies is that their 
stellar popul ations are relativ ely simple, which can possibly be modeled by two population 
components (ILi et al.l l2012a] ). Our main purpose is to study how the results of spectral 
fitting can be changed by using bsSP models instead of ssSP models. UV-optical spectra 
will be fitted by homogeneous ssSP and bsSP models that vary from only including binary 
interactions or not when evolving stars. 

The structure of this paper is as follows. Sect. 2 introduces stellar population models, 
and Sect. 3 describes spectral synthesis. Sect. 4 then presents the spectral fitting code 
(BS2fit) used in this work. In Sect. 5, for some mock galaxies, the galaxy parameters 
determined by ssSP and bsSP models are compared. Then in Sect. 6, a sample of 10 early- 
type galaxies are used to give a similar test. Finally, Sect. 7 summarizes and discusses on 
the results. 



2. Stellar population model 



A n advanced stell ar population synthesis model including both ssSPs and bsSPs flLi fc Han 



2008al : iLi et al.ll2012al ) is taken for this work. This model supplies homogeneous spectra for 



both simple and composite stellar populations (SSPs and CSPs) with or without binary 
interactions. This makes it possible to study the result changes caused only by binary evo- 
lution. In the following two subsections, we introduce the stellar population synthesis model 
for SSPs and CSPs, respectively. 



2.1. Simple stellar populations 



Simple stellar populations (SSPs) are built as follows. Firstly, some sta rs (primary com- 
ponents of binaries) are generated following an initial mass function (IMF) of lChabrierl (j2003[ ) 
with lower and upper mass limits of 0.1 and 100 M0 respectively. Then for each binary, the 



-4- 



mass of its secondary component is generated by taki ng an uniform distribution for the mass 
ratio (g, 0-1) of secondary to primary component (IMazeh et al.l Il992l : iGoldberg fc Mazeh 
1994I ). The separation (a) between two binary components is given following an assumption 
that the fraction of binary in an interval of log(a) is constant whe n a is big (lO-Rry^ < a < 
5.75 X lO^i?©) and it falls off smoothly when a is small (< lOi?©) (IHan et al.lll995l ). which 
can be written as 



a .p{a) 



"Sep 



[a/ao) 



where a. 



Sep 



0.070, ao = lOi?0,ai 



isep, aQ < a < ai 

: 5.75 X lO^i?© and i) ^ 



(1) 

1.2. The eccentricity (e) 
of each binary system is generated according to an uniform distribution. In this method, 
each population contains about 50% binaries with orbital period less than 100 yr, which 
is similar to the typical binary fraction of the Galaxy. When building binary star simple 
stellar populations (bsSSPs), the interactions between two binary components are taken into 
account when evolving stars, but all stars are evolved separately when building single star 
simple stellar populations (ssSSPs). All stars of an SSP are assumed to form in a star burst 
and have the same metallicity. 

After the generation of po pulation stars , all o f them (4 000 000) are evolved using the 



rapid stellar evolution code of iHurley et al.l (|2002[ ) (Hurley code). Most binary evolution 
processes such as mass transfer, mass accretion, common-envelope evolution, collisions, su- 
pernova kicks and angular momentum loss are included for bsSPs. Different mass trans- 
fer mechanisms, i.e., dynamical mass transfer, nuclear mass transfer a nd thermal mass 
transfer are taken into acc ount using the results of rnany w orks (e.g., iTout et al.l Il997l : 
Hjellming fc Webbinklll987l ). One can see iHurley et al.l ( 12002| ) for more details. Some de- 
fault values for Hurley code, i.e., 0.5, 1.5, 1.0, 0.0, 0.001, 3.0, 190.0, 0.5, and 0.5, are taken for 
wind velocity factor (/?„), Bondi-Hoyle wind accretion faction (ctw), wind accretion efficiency 
factor (/iw), binary enhanced mass loss parameter (B^), fraction of accreted material retained 
in supernova eruption (e), common-envelope efficiency (acE)? dispersion in the Maxwellian 
distribution for the supernovae kick speed (uk), Reimers coefficient for mass loss (77), and 
binding energy factor (A), respective l y. We take these default values, because they have been 
checked in the work of iHurley et al.l (120021 ). Although these default values remain somewhat 
large uncertainties, the results for spectral stellar population synthesis will be not affected 
too much by the uncertainties in these parameters, according to a test in our previous work 
( jLi et al.ll2012al ). In addition, because Hurley code uses some fitting formulae to calculate 
the evolution of stars, it causes about 5% uncertainty in the evolutionary parameters of stars. 

The evolutionary parameters of stars are finally transformed to the spectra l energy distri- 



butions (SEP) (or spectra) of ste llar populations by BaSeL 3.1 spectral library ( iLejeune et al. 



1997, 1998; Westera et al 



20021). The library is chosen here because of its wide wavelength 
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coverage and reliability. The uncertainties in the spectra of populations caused by the spec- 
tral library is small (about 3% on average). In the same way, the homogenous spectra of 
both ssSSPs and bsSSPs are computed. Note that the only difference between two kinds of 
models is that bsSSPs take binary interactions into account but ssSSPs do not. 

When comparing the spectra of bsSSPs to those of ssSSPs, we find obvious difference. 
The left panels of Fig. 1 show the comparison of spectra of two kinds of populations. It 
is clear that the spectra of bsSSPs and ssSSPs are different, especially in UV band. Many 
old (> 3Gyr, bottom lines) bsSSPs show U V-upturn spectra , which is mainly caused by 



hot subdwarf and blue straggler stars (e.g., iHan et al.l 120071 ). but all old ssSSPs do not 



have sim ilar spectra. Bec ause UV-upturn phenomenon has been observed in many elliptical 
galaxies ( IO'Connelllll999l ). which are usually thought as SSPs, our result suggests that bsSSPs 
can better fit to early-type galaxies than ssSSPs. In binary evolution, the formation channels 
for hot subdwarfs are common-envelo pe ejection, stab le Roche lobe overflow, and merger of 
helium white dwarfs (WDs) (see also lHan et al.ll2007l ). 



2.2. Composite stellar populations 



The SEDs of CSPs are built on the basis of SEDs of SSPs. Because there is no common 
result for the star formation histories (SFHs) of early-type galaxies, a simple method is 
taken to model CSPs. In detail, each CSP is assumed to contain a pair of old and young 
components with the same metallicity. This assumption is in agreement of previous studies 
that early-type galaxies are dominated by old populations and there is only a little fraction 
of young populations in such galaxies. The mass fraction of young component is assumed to 
be dependent on the ages of two components of CSPs, which is calculated by formula (2). 
It means that the mass fraction of young component declines exponentially with increasing 
difference between the ages of old and y oung componeii t s. Th is agrees with previous studies 



on the SFHs of early- type galaxies, e.g., iThomas et al.l ( 120051 ) 



F2 = 0.5 exp[ 



to-U 



(2) 



where F2 is the mass fraction of young component; ti and ^2 are the ages of old and young 
components of a CSP, respectively. As a standard model, r is taken as 3.02, according to 
the observational fraction o f bright early-type galaxies wit h recent (< 1 Gyr) star formation 
at a level more than 1-2% JyI et al.lboosi : Ihi fc Hani [20071 ). 



By the above method, the SEDs of single and binary star composite stellar populations 
(ssCSPs and bsCSPs) are calculated. In the right panels of Fig. 1, the SEDs of two kinds 
of CSPs are compared. For some pairs of ssCSP and bsCSP with the same parameters. 
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we see clear difference between the SEDs of two kinds of populations. Similar to the case 
of SSPs, the UV spectra of two kinds of populations are obviously different. In addition, 
we find that both ssCSPs and bsCSPs can show UV-upturn spectra, but UV-upturn comes 
from different reasons. The UV-upturn of bsCSPs are mainly caused by binary evolution, 
but that of ssCSPs by young stars. This suggests that bsCSP and ssCSP models possibly 
give different estimates for the SFHs of early-type galaxies. 



3. Spectral synthesis 



The observed spectra of stellar populations are built based on the SEDs of CSPs. 
According to our model, each CSP is described by three parameters: metallicity Z, old- 
component age ti, and young- component age t2- The effects of stellar velocity dispersion 
cr*, extinction from dust around stellar population, redshift, and the extinction caused by 
the Milky Way are added to the SEDs of CSPs to form the observed spectra. Similar to 
the works of ICid Fernandes et al.l (120051 ). the Line-of-sight stellar motions are modeled by a 
Gaussian distribution G centered at velocity t>* and with dispersion cr*. The effect of dust 
around pop ulation, which is param etrized by y-band optical depth ry, is modeled following 
the work of ICharlot &: Fall! (120001 ). When modeling the Gal actic extinction, whi ch can be 
parametrized by \^-band extinction Ay, the extinction law of ICardelli et al.l (Il989[ ) with Ry 
= 3.1 is adopted. The observed flux (fx) at wavelength A can be expressed by 



fx = [fxo{z,ti,t2) ® G(w*,a^)] p{tx) r{Ax), 



(3) 



where fxo is the flux of a CSP with three parameters {Z,ti,t2), ^(f*, a^:) is a Gaussian dis- 
tribution with mean and standard deviation of f * and a* respectively, p{tx) is the percentage 
of energy that pass through the dust with optical depth of tx, r{Ax) is the fraction of energy 
after the extinction (Ax) of Milky Way galaxy. ® denotes the convolution operator. 

Fig. 2 gives two examples for the synthesis of the spectra of a pair of ssCSP and bsCSP. 
We see that Line-of-sight stellar motions (redshift and stellar velocity dispersion) affect the 
wavelength range and strength of absorption lines obviously, but stellar velocity dispersion 
does not affect continuum spectra too much. Meanwhile, the dust around stellar population 
and the Milky Way affects the whole spectra. We can also see significant difference between 
the spectra of ssCSPs and bsCSPs from this figure. 
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Fig. 1. — Comparison of SEDs of single and binary star populations. Grey lines are for 
single star populations, while red lines for binary star populations. Left and right columns 
show the results of simple and composite stellar populations (SSPs and CSPs), respectively. 
Note that stellar age increases from top to bottom in each panel. 
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Fig. 2. — Formation of synthetic SEDs of single (left) and binary (right) star composite stellar 
populations (ssCSP and bsCSP). Green lines show SEDs of two simple stellar populations 
(SSPs) with age = 10 Gyr and solar metallicity {Z = 0.02), while black, red, blue, cyan, 
and purple lines are SEDs after adding the effects of second star formation, stellar velocity 
dispersion, dust in galaxy, redshift, and Galactic extinction into the original SEDs of two 
SSPs, gradually. 
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4. Spectral fitting: BS2fit code 



Spectral fitting is an important technique to determine many properties of galaxies, 
e.g., redshift, stellar metallicity, age, velocity dispersion, mass, and dust extinction, because 
different wavelengths are usually dominated by various physical processes. Thus it is a good 
choice to study the role of binary evolution in stellar population studi es via spectral fitting . 



Although th ere are some available spectral fitting c odes, e.g., MOPED (Heavens et al.ll2000f) 



(ITojeiro et alj 



2007 



PLATEFIT (ITremqnti et al.ll2004[) . STAR LIGHT (ICid Fernandes et al 



20091). STECKM AP flOcvirk et all 120061). sedfit (IWalcher et all 120061 ) 



2004l2005h. V ESPA 



NBURSTS (IChilingarian et all 120071 ). and ULYSS flKoleva et all 120091 ). they are not very 
suitable for this work. Firstly, all of them are developed for some special works so that 
their fitting methods and wavelength coverages are not suitable for this work. For example, 
STARLIGHT aims to fit the optical spectra of SDSS galaxies but our work needs a wider 
wavelength coverage. VESPA studies the SFHs of galaxies via some binned parameters, but 
we want to do more accurate studies. Furthermore, there are some other limitations to use 
these codes, e.g., the dependence of internet connection or unaltered procedures. Therefore, 
we decide to create a new code. We call our new code binary star to fit (hereafter BS2fit), 
which can be used for spectral fitting via both bsSPs and ssSPs, within a wide wavelength 
range from UV to optical band. The new code can be briefly introduced as follows. 

BS2fit code will be revised and updatecll]. At this stage, BS2fit aims to determine a 
few parameters of early-type galaxies from observed spectra using bsSP or ssSP models. An 
observed spectrum is assumed to be determined by eight input parameters, i.e., Z, ti, t2, f*, 
cr*, Tv, Ay, and total stellar mass M*. Because there are so many input parameters, it is 
not easy to get a best-fit result by directly searching for the combination of parameters with 
minimum in a high-resolution parameter grid. Therefore, we divide the fitting task into 
three parts. First, Line-of-sight velocity (or redshift) is determined by a few methods that 
are widely used (e.g., a bsorption or emission line fitting and cross-correlation analysis) (see 
e.g.. ICoUess et al.ll200ll ). Then the Galactic extinctio n parameter Ay is determined from the 



direc tion of galaxy and dust map of the Galaxy (e.g., iBurstein fc Heileslll982[ ISchlegel et al. 



19981 ). Finally, the other parameters are determined based on fixed redshift and Galactic 
extinction. Although this has made fitting process much faster, it is still difficult to get 
the result by a direct search on a personal computer. Two techniques are therefore used to 
make the fitting significantly quick. One is to estimate the possible ranges of parameters 
using SSP models before comparing the observed spectrum to those of CSPs. We call this 



^The code and data can be obtained on request to Zhongmu Li, and the authors are trying to make it 
available on the internet. 



- 10 - 



procedure SSP-fitting. This is very effective, because it guides the code to search in only 
hmited ranges of parameters. This method is also reliable for most populations, because SSPs 
usually estimate lower metallicities, smaller stellar velocity dispersions, and younger ages for 
populations, compared to CSPs. In detail. Fig. 3 shows how the ranges of parameters of 
74 test stellar populations are related to SSP-fitted parameters. On average, the parameter 
ranges of about 95 percent of CSPs can be well estimated from SSP-fitting. 

The other technique for speeding up the fitting is to compare continuum spectra before 
comparing whole spectra. This skips many populations according to their continuum differ- 
ence from the observed spectra. This can be done because there is a maximum change in 
continuum spectra corresponding to the effects of dust extinction and stellar velocity disper- 
sion when their ranges are given. The goodness of a fit is judged by = ^[{fohx — fthx)"^ ^x], 
where fohx and fthx are observed and theoretical spectra, while u\ is the weight for wave- 
length A. The best-fit result corresponds to minimum x^- 

We mask out some emission li ne regions as they can affect the final results obviously. 
Acco r ding to some previous works (IMcQuade et al.l Il995l : iPanter et al.l 120071 : iTojeiro et al. 
20071 . 120091 ). the following emission- line regions in every spectrum's rest frame wavelength 
range are masked out: 1213-1219, 1548-1550, 1637-1643, 1905-1911, 2323-2329, 2421-2427, 
2796-2803, 3423-3429, 3711-3741, 4087-4117, 4325-4355, 4846-4876, 4944-4974, 4992-5022, 
5885-5900, 6535-6565, 6548-6578, 6569-6599, 6702-6732, 6716-6746, 6728-6734, 7132-7138, 
7319-7330 A. 

Because we want to make the results accurate enough, a high-resolution grid of parame- 
ters is taken. In detail, stellar metallicity (Z) is from 0.0003 to 0.03, with an interval of 0.0001 
when Z is less than 0.001 and an interval of 0.001 for higher metallicity. The old-component 
age (ti) of populations is in the range of zero to 15 Gyr, with an interval of 0.1 Gyr, and the 
young- component age (^2) decreases from ti to zero. Stellar velocity dispersion, cr*, changes 
from to 350 km s~^ by a step of 10km s~^, and the optical depth (tv) is between and 
1.5, with an interval of 0.01. Note that the range s of and Ty are chosen according to some 
previous works, e.g., ICid Fernandes et al.l ( 120051 ) and lTojeiro et al.l (|2007| ). 



In order to check the accuracy of BS2fit code. Figs. 4 and 5 show the input versus 
recovered values of mass-weighted age < >m, luminosity-weighted age < >l, Z, ti, 
t2, "W*, 0"*, Ty, and M^: of 37 ssCSPs (open squares) and 37 bsCSPs (filled squares). In the 
test relating to Fig. 4 all wavelengths are assumed the same weight because there is no 
error in spectra, while in the study corresponding to Fig. 5 the weight of each wavelength 
is set to be l/E"^-)^, where Eox is observational error. The stellar mass of each population is 
calculated from the input flux at 5500 A and mass-to-light ratio of best-fit population. We 
see that when there is no error in input spectra, BS2fit code can well recover most parameters 
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Fig. 3. — Comparison of SSP-fitted and real parameters of 37 ssCSPs and 37 bsCSPs. Open 
and filled squares are for ssCSPs and bsCSPs, respectively. Solid lines show the ranges for 
searching the final parameters via CSPs, while dashed lines are unity lines, ti is in Gyr and 
(7* is in km s~^. 
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(Fig. 4), although the t2 values of some populations with large (> 3.5 Gyr) input t2 are not 
recovered accurately because of the degeneracy effects of various parameters on the same 
spectra. Note that the input spectra of a pair of ssCSP and bsCSP with the same set of 
input parameters are usually different. 

If there are uncertainties (signal-to-noise ratio S/N = 10) in input spectra (Fig. 5), 
BS2fit code reports obviously larger ages (especially ti) for most populations. This because 
the change of spectra is less sensitive to the age change of old populations. Moreover, the 
uncertainties in input spectra result in lower estimates for stellar metallicities of some pop- 
ulations, because the change of spectra has different sensitivities to the metallicity changes 
of metal-poor and rich populations. The deviations in stellar age and metallicity are also 
related to the well-known age-metallicity degeneracy. Therefore, the S/N of observed spec- 
tra is crucial for constraining the stellar metallicities and ages of galaxies. In addition, we 
see that usually bsCSP and ssCSP models lead to similar difference between the input and 
recovered parameters. This means that BS2fit code can be used to study the difference 
between galaxy parameters, which are caused by taking various (bsCSP and ssCSP) models, 
although the values of parameters may be not so accurate. Note that the fitted parameters 
and uncertainties are given statistically by taking the average and 1 a deviation of the results 
of fitting some pertur bed spectra that are generated randomly accor ding to the observational 



uncertainty (see also lCid Fernandes et al.ll2005l : iTojeiro et al.ll2007l ) 



5. Results from mock galaxies 

This section studies the difference between bsSP and ssSP-fitted results, via a sample of 
37 mock galaxies, which covers wide parameter ranges. We assume that all mock galaxies are 
bsCSPs, as almost all galaxies may contain some binaries. In addition, because any observed 
spectrum possibly contains some uncertainties, the effect of observational uncertainties is 
taken into account. Similar to the previous section, the S/N of spectra is set to 10, and u\ 
is given as 

Fig. 6 shows the comparison of bsSP and ssSP-fitted parameters of 37 bsCSPs (mock 
galaxies). We see that ssSP and bsSP- fitted results for mass and luminosity- weighted ages, 
metallicity, stellar mass, old-component age, stellar velocity dispersion, and galactic dust 
extinction are roughly consistent when uncertainties are taken into account. In particu- 
lar, ssSP and bsSP-fitted stellar masses are almost the same. However, bsSP-fitted young- 
component ages (t2) are different from ssSP-fitted results, significantly. For mock galaxies 
with ssSP-fitted t2 younger than 3.5 Gyr, bsSP models measure much older young popula- 
tion components for most mock galaxies. This suggests that ssSP models usually measure 
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Fig. 4. — Comparison of input and recovered parameters of ssCSPs and bsCSPs. The input 
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more recent star formations for early-type galaxies, compared to bsSP models. In this case, 
m aybe many early-ty pe galaxies do not contain so much recent star formations as the results 
of I Chen et al.l (120 12[ ). which was obtained using some ssSP models. 



6. Results from real galaxies 



In this section, we try to compare the galaxy parameters that are determined via ssSP 
and bsSP models, for 10 early-type galaxies. The data are selected from a catalogue of 
UV-optical spectra of 99 local galaxies, which is available on the internet!^. The advantage of 
using these data is the wide wavelength coverage from UV to op t ical band. Beca use binaries 
affect UV spectra obviously (IHan et al.l 120071 : iLi fc Hanll2008al : iLi et al.ll2012al ). such data 
seems ideal for exploring the effect of binary evolution on spectral fitting. Besides some 
emission lines, we also mask out the spectra around two conjunction wavelengths (2000 and 
3000 A) as the two parts contain obviously larger uncertainties. In addition, because the 
quality of optical spectra is better than UV spectra, optical spectra are given higher weights 
by taking S/N = 8, while the S/N of UV spectra is taken as 7. The S/N values are given 
according to a description of the data (IKinney et al.lll996l ). but the data quality seems not 
as good as S/N = 8, because the minimum values of fittings are obviously larger than 1. 

Using BS2fit code, six parameters of 10 galaxies are determined, and the ssSP and 
bsSP-fitted results are then compared. For convenience, we take the redshifts of 10 galax- 
ies ftwi_^£ASA7lPA^ Extrag^^ (NED). The origi n al results were obtaine d 
bvlWegner et al.l feoosh: Ida Costa et al.l fll99lh : I tuIIv et al.l ^A: korib alski et al.l f|20oi): 



Huchra et al. 



( 119991 ): iDenicolo et al.l (120051 ). The Galactic extinction s are calculated usin g 
the extinction calculator of NED, which uses the data and technique of ISchlegel et al.l ( 1l998l ). 
Tables 1 and 2 list the ssSP and bsSP-fitted parameters of 10 galaxies. Because the qual- 
ity of observed spectra is not very good, the result for stellar velocity dispersion may be 
not accurate. Fig. 7 compares the best-fit spectra of ssSP and bsSP fittings with observed 
ones, and Fig. 8 compares ssSP and bsSP-fitted parameters. We find that both ssSP and 
bsSP models can fit to the spectra of eight galaxies, i.e, NGC1399, NGC1553, NGC2865, 
NGC3031, NGC1404, NGC221, NGC1052, and NGC205. At the same time, both bsSP and 
ssSP models do not fit the spectra part with wavelength larger than 6000 A for NGC210. 
This may result from the low upper metallicity (0.03) of our theoretical populations. In 
addition, we find that ssSPs can fit the observed spectra of NGC1433 better than bsSPs. 
This implies that most stars of this galaxy may be single stars, rather than binary stars. 



^ftp:/ /ftp. stsci.edu/pub/catalogs/nearby_gal/sed. html 
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Fig. 6. — Comparison of ssSP and bsSP-fitted parameters of 37 bsCSPs. Points and error 
bars show mean parameters and a uncertainties, respectively, "ss-fitted" and '"bs-fitted" 
denote ssSP-fitted and bsSP-fitted results, respectively. 
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Fig. 8 shows the comparison of ssSP and bsSP-fitted parameters of 10 galaxies. We see 
that ssSP and bsSP models report obviously different results for many parameters of galaxies. 
In particular, the bsSP-fitted young-component ages (^2) of six galaxies are significantly 
older than ssSP-fitted results, and the bsSP-fitted dust extinctions (r„) of half galaxies are 
obviously larger than ssSP-fitted results. In addition, two kinds of stellar population models 
give different results for the mass-weighted, luminosity-weighted, and old-component ages 
of some galaxies. Furthermore, we see some difference between the bsSP and ssSP-fitted 
stellar velocity dispersions (a*). However, bsSP and ssSP models measure similar stellar 
metallicities and masses for most galaxies. When comparing the results of 10 real galaxies and 
37 mock galaxies, we find good consistency: bsSP models measure older young-component 
ages and larger dust extinctions than ssSP models, while two types of models measure similar 
stellar metallicities and masses for galaxies. Because BS2fit can recover the ages of recent 
(^2 < 3.5 Gyr) star bursts well (see Fig. 4), both tests of mock and real galaxies indicate 
that binary evolution affects the measurement of SFHs of early-type galaxies, significantly. 



7. Conclusion and discussion 

This paper investigates the potential importance of binary evolution in UV-optical spec- 
tral fitting of early- type galaxies, via a new spectral fitting code (BS2fit). Both the results 
of mock and real galaxies show that stellar population models including binary evolution 
will give different estimates for galaxy parameters such as young- component age and dust 
extinction. The difference between bsSP and ssSP-fitted young-component ages {t2) can be 
as large as 7 Gyr (mock galaxies) or 10 Gyr (real galaxies), while the optical depth (ry) can 
be different by about 0.15 if ssSP-fitted ry is low (< 0.25). Our results suggest that binary 
evolution can possibly play an important role in spectral fitting of early-type galaxies. It is 
therefore necessary to take binaries into account in spectral synthesis studies. Because stars 
of most galaxies are possibly formed in more than one star bursts, and SFHs supply unique 
information for studying galaxy formation and evolution, it is crucial to include the effects 
of binaries in stellar population models when studying the SFHs of galaxies. In addition, 
this study shows that bsSP models measure larger optical depth for galactic dust, because 
bsSPs are more luminous than ssSPs in UV band. This is possibly helpful for explaining why 
some minus extinction valu es are needed in some spectral fittings based on ssSP models (e.g.. 



Cid Fernandes et al.l 120051 ). Furthermore, although parameters such as old-component age 
and stellar velocity dispersion, which are derived from bsSP and ssSP models, are roughly 
consistent, the results for some specific galaxies are different a lot (e.g., 8 Gyr for ti and 100 
km s~^ for o"*). This implies that binary evolution may affect more studies about galaxies. 
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Fig. 7. — Comparison of bsSP and ssSP-fitted SEDs with observational spectra of 10 galaxies. 
Black, red and blue lines are for observed, bsSP-fitted, and ssSP-fitted spectra, respectively. 
All spectra are normalized at 5500 A. 
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Fig. 8. — Comparison of ssSP and bsSP-fitted parameters of 10 galaxies. Solid lines are 
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respectively. 
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Table 1. ssSP-fitted parameters of 10 galaxies. The unit oiti,t2, < >m and < >l is 
Gyr, and that of cr* is km s""^. The values of redshift and Ay are taken from NED. 



Galaxy 


redshift 




Ay 




z 






tl 




t2 








TV 




< t* >M 


< t* >L 


NGC1399 


0.004753 


0, 


.0393 


0, 


.0263±0. 


.0015 


13. 


.55±0. 


.14 


4.92±6. 


.59 


237±195 


0. 


.138±0. 


.097 


13.38±0. 


.24 


12 


.63±0.61 


NGC1553 


0.004236 


0, 


.0468 





.0172±0. 


.0035 


9. 


.75±1. 


.65 


2.40±3, 


.13 


333±44 


0, 


.143±0. 


.172 


9.46±1. 


.76 


8 


.42±1.49 


NGC2865 


0.008619 


0, 


.1608 


0, 


.0172±0. 


,0040 


10. 


18±3. 


,92 


1.30±0, 


.85 


350±0 


0, 


.037±0. 


,050 


9.90±4. 


15 


8 


.62±3.75 


NGC3031 


-0.000113 


0, 


.1486 


0, 


.0257±0. 


,0006 


13. 


,50±0. 


,00 


0.82±0, 


.05 


350±0 


0, 


.382±0. 


,029 


13.40±0. 


,00 


12 


.19±0.06 


NGC1404 


0.006391 


0, 


.0353 


0, 


.0205±0. 


,0036 


13. 


,45±0. 


,25 


0.80±0, 


.16 


314±68 


0, 


.130±0. 


128 


13.35±0. 


,26 


12 


.26±0.11 


NGC 210 


0.005457 


0, 


.0673 





.0160it0. 


.0061 


6. 


.68it0. 


.72 


1.87±0. 


.23 


342±22 


0. 


.146it0. 


.183 


6.20±0. 


.75 


5, 


.28±0.80 


NGC 221 


-0.000667 


0, 


.4797 





.0113±0. 


.0074 


11. 


.25±4. 


.04 


1.38±0. 


.20 


350±0 


0. 


.124±0. 


.197 


11.02±4. 


.22 


10 


.00±4.57 


NGC1433 


0.003586 


0, 


.0282 





.0155±0. 


.0052 


6. 


.20±0. 


.00 


3.78±2, 


.74 


350±0 


0, 


.102±0. 


.120 


5.91±0. 


.26 


5 


.46±0.83 


NGC1052 


0.001004 


0, 


.0824 


0, 


.0200±0. 


.0000 


13. 


.60±0. 


.00 


0.58±0, 


.05 


350±0 


0. 


.138±0. 


.024 


13.51±0. 


.00 


11, 


.79±0.15 


NGC 205 


-0.000777 


0, 


.1685 





.0122±0. 


.0005 


6. 


.20±0. 


.00 


1.60±0. 


.00 


307±64 


0. 


.OOOitO. 


.000 


5.70±0. 


.00 


5 


.02±0.01 



Table 2. Similar to Table 1, but for bsSP-fitted parameters. 



Galaxy 


redshift 


Av 


Z 


h 


t2 




TV 


< U >M 


< >L 


NGC1399 


0.004753 


0.0393 


0.0300±0.0000 


13.50±0.00 


12.27±0.81 


350±0 


0.182±0.062 


13.13±0.20 


13.10±0.23 


NGC1553 


0.004236 


0.0468 


0.0183±0.0069 


11.87±2.56 


1.18±0.19 


348±5 


0.265±0.090 


11.69±2.67 


10.55±2.93 


NGC2865 


0.008619 


0.1608 


0.0145±0.0033 


14.48±0.30 


0.68±0.16 


348±5 


0.043±0.060 


14.41±0.31 


12.98±0.29 


NGC3031 


-0.000113 


0.1486 


0.0230±0.0000 


12.80±0.00 


10.53±0.97 


350±0 


0.337±0.019 


12.29±0.04 


12.20±0.11 


NGC1404 


0.006391 


0.0353 


0.0197±0.0045 


11.33±0.45 


3.25±3.43 


348±5 


0.232±0.075 


11.07±0.52 


10.33±0.56 


NGC 210 


0.005457 


0.0673 


0.0173±0.0006 


6.20±0.00 


5.95±0.14 


350±0 


0.235±0.115 


6.09±0.06 


6.09±0.06 


NGC 221 


-0.000667 


0.4797 


O.OllOiO.OOOO 


7.60±0.00 


4.43±2.19 


350±0 


0.165±0.058 


7.18±0.19 


6.85±0.44 


NGC1433 


0.003586 


0.0282 


0.0172±0.0043 


6.67±1.27 


3.68±3.01 


350±0 


0.070±0.097 


6.44±1.28 


5.66±0.71 


NGC1052 


0.001004 


0.0824 


0.0200±0.0000 


ll.lOiO.OO 


ll.lOiO.OO 


350±0 


O.OOOiO.OOO 


ll.lOiO.OO 


ll.lOztO.OO 


NGC 205 


-0.000777 


0.1685 


0.0032±0.0005 


14.70±0.49 


0.20±0.00 


321±51 


0.013±0.036 


14.64±0.50 


9.95±0.69 



A fixed binary fraction (50%) is taken for all bsSPs, although the real fraction may 
be different from a galaxy to another. However, our results are potentially useful, because 
50% is a typical binary fraction of our Galaxy, and this fraction can fit the spectra of most 
galaxies. However, the assumption of separation of binary components may have obvious 
effects on the final results, althoug h other adjustable parameters affect the results slightly 
(IHurley et al.ll2002l : iLi et al.ll2012al ). On another side, although it is possible that UV- upturn 
is caused mainly by binary evolution, some single ext reme horizontal br anch (EHB) stars 
can also affect the UV flux of early-type galaxies (e.g., iBrown et al.l 120001 ). Thus our results 
actually show the maximum effects of binaries on spectral fittings of early-type galaxies. In 
addition, a simple assumption for SFH was used to model the CSPs of early-type galaxies. 
This is possibly reasonable for most galaxies, but it is interesting to do some researches via 
taking different SFHs. Furthermore, the stellar masses of populations were measured by the 
flux at 5500 A in this work. If stellar mass is determined at other wavelengths, e.g., 4000 A, 
the bsSP and ssSP-fitted stellar masses may be different more. Moreover, we used a rapid 
spectral fitting technique in BS2fit code because of the limitation of our computation ability. 
This leads to some uncertainties in the results. Although this does not affect the results of 
this work too much, it is necessary to make the code more accurate and check the results 
again. 
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